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ABSTRACT 



Recent observations and modeling of globular clusters with multiple populations 
strongly indicate the presence of super helium-rich subpopulations in old stellar sys- 
tems. Motivated by this, we have constructed new population synthesis models with 
and without helium-enhanced subpopulations to investigate their impact on the UV- 
upturn phenomenon of quiescent early- type galaxies. We find that our models with 
helium-enhanced subpopulations can naturally reproduce the strong UV-upturns ob- 
served in giant elliptical galaxies assuming an age similar to that of old globular clusters 
in the Milky Way. The major source of far-UV (FUV) flux, in this model, is relatively 
metal- poor and heliu m-enhanced hot horizontal branch stars and their progeny. The 
Burstein et al.l (|1988l ) relation of the FUV — V color with metallicity is also explained 
either by the variation of the fraction of helium-enhanced subpopulations or by the 
spread in mean age of stellar populations in early-type galaxies. 



Subject headings: galaxies: elliptical and lenticular, cD - galaxies: evolution — galaxies: 
stellar content - ultraviolet: galaxies — galaxies: individual (M87) 



1. INTRODUCTION 



It is well established that the far ultraviolet (FUV) flux ( "UV-upturn" ) observed in the nearby 
quiescent early-type galaxies (ETGs) origi nates from a minority population of hot hori z ontal-branch 



(HB) stars and their progeny (see, e.g., iGreggio Renzinil I199Q : lO'Connelll Il999i : iBrown et al 
2000). Recently, a significant progress has been made as to the origin of these hot HB stars 
in old stellar systems. This new theory is based on the recent observations and modeling of 
globular clusters (GCs) with extended HB, such as u Cen and NGC 2808, where the multiple main- 
sequences and h ot HBs in these clusters can only be explained by the presence of super h e lium- rich 



subpopulations (jD'Antona k, Caloi 



2004 



Norris 



2004 



Lee et al 



2005 



Piotto et al. 



2005 



20071 ). 



The origin of this helium enhancement is most likely due to the pollution from the intermediate- 
mass asymp totic giant branch stars and/or fast-rotating massiye stars, or due to the enrichment b y 



super novae (j Ventura &: D'Antona 



2008 



Decressin et al. 



2007 



Piotto et al. 



2005 



Lee et al J 12003 ). 
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Therefore, the likely presence of super-helium-rich subpopulations and the resulting hot HB stars 
in old stellar systems deserve further investigation, as they could be a major source of the FUV 
flux in quiescent elliptical galaxies. The purpose of this Letter is to report our first result on the 
UV upturn phenomenon predicted from the Yonsei evolutionary population synthesis models with 
helium-enhanced subpopulations. 



POPULATION SYNTHESIS MODELS 



The models presented in this paper were construct ed using the Yonsei E volutionar y Popu lation 
Synthesis (YEPS) code. The readers are referred to iPark Led (jl997l ). iLee et al.l (|2000l ). and 
Chung et al. (2011, in prep.) for the details of model construction. In order to investigate the 
effect of helium-enhanced stellar population on the UV-upturn, we have used the most up-to-date 
Yonsei- Yale (Y 2 ) stellar isochrones and HB evolutionary tracks (Han et al. 2011, in prep. ) with 
different values of helium abundance (Y = 0.23, 0.33, and 0.38). We apply the same lReimerd (|1977l ) 



mass- loss parameter (r/) and the enhancement in a-elements ([a/Fe]) for both helium-enhanced 
and normal helium populations (see Table 1). 

Figure 1 shows an example of our model HR diagrams and corresponding spectral energy 
distributions (SEDs) of simple stellar population (SSP) for normal helium and helium-enhanced 
populations. Since helium-rich stars evolve faster than helium-poor stars, helium-rich stars have 
lower masses at given age. This effect has a striking difference on the HB stage (see left panels). The 
mean temperature of HB stars in helium-enhanced (Y = 0.33) population is ~ 11,500 K greater 
than that of the normal helium (Y = 0.23) population at a given metallicity ([Fe/H] = —0.9) 
and age (11 Gyr). Consequently, the helium-enhanced models show extremely strong far-UV flux 
compared to the normal helium models. 

In order to include this effect of the helium-enhanced population into the composite SED 
for a model elliptical galaxy, we have employed following procedures in the model construction. 
First, helium abundances of helium-enhanced subpopulations were chosen to be Y = 0.33 and 
0.38, as these values are roughly required to reproduce the extreme blue HB stars observed in 
the Milky Way GCs (|p'Antona fc Caloilbooi boosl : Lee et al.l I2OO5I : IPiotto et al.1 1200.4 Joo & Lee 
2011, in prep.). In our modeling, populations having these two helium abundances were mixed 
half and half. Second, we have used mean FUV — V colors of UV-strong GCs in M87 (jSohn et al 



2006) to set the fraction of helium-enhanced subpopulations at given metallicity (see Figure 2a). 
Figure 2b shows the fraction of helium-enhanced subpopulations as a function of metallicity, which 
was adopted in our model with helium-enhanced subpopulations to reproduce the observed trend 
of FUV — V color with metallicity in M87 GCs. Figure 2c displays the contribution of helium- 
enhanced subpopulations in the metallicity distribution function (MDF) used in our construction 
of the composite stellar population mod el for giant elliptical g alaxie s. This MDF was adopted from 
the simple chemical evolution model of Kodama fc Arimoto (jl997l ). 
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According to the list of lLee et al.1 (|2007l ). about 30% of the Milky Way GCs have extended HB 



each of wh i ch contains about 30% of helium- e nhanced subpopul ation (|Lee et al.ll2005l ; iPiotto et al 



20051 . l2007l : lYoon et al.ll2008l : lHan et al.ll2009l : Bellini et al.ll2010l ). so about 9% of stellar population 
in the Milky Way GCs are assumed to be helium-enhanced population. In our composite model for 
giant elliptical galaxies, the number fraction of helium-enhanced subpopulations is about ~ 11%, 
which agrees roughly with that estimated in the Milky Way GCs. This fraction of helium-enhanced 
subpopulations is reduced to ~ 6.4% when only hot enough (T e ff > 20, 000 K) blue HB stars are 
counted, which are actually responsible for the far-UV flux iaour model giant elliptical galaxiesf] 



COMPARISON WITH OBSERVATIONS 



Figure 3 presents our composite models constructed with and without helium-enhanced sub- 
populations compared with the SEDs observed by IUE satellite for two giant elliptical galaxies 
NGC 4552 and NGC 4649. In our models, the value of mean metallicity ( ( [Fe/H] ) ) was chosen so 
that they can reproduce the observed Mg b 4.6 A) i ndex of theses galaxies (see also Figure 4) , 
assuming that [a/Fel is 0.3 in giant elliptical galaxies (jWorthey et al.lll992l : iThomas et al.ll2005l : 
Kormendv et alihood ). The parameters adopted in our models presented in Figure 3 are listed in 
Table 1. It is clear from Figure 3 that our models with helium-enhanced subpopulations can natu- 
rally reproduce not only the strong far-UV upturns but also the 2500 A dips observed in NGC 4552 
and NGC 4649, while the models without helium-enhanced subpopulations fail to repro duce the far- 
UV upturns, unless the age of underlying stellar population is increased to ~ 16 Gyrs (IPark k, Lee 
19971 ; lYi et al.lll999l ). Therefore, the minority population (~ 6.4%) of helium-enhanced hot HB 
stars and their progeny is responsible for the observed far-UV flux in this scenario!! 

Our composite models with helium-enhanced subpopulations predict that about 85% of the 
far-UV flux comes from the metal-poor side ([Fe/H] < 0. 0) of the MDF. In this regard, this model 
is qualitatively similar to the "metal-poor HB model" of Park &: Lee! (1997), but it does not need 
to inv oke unrealistically old ages (~ 16 Gyr) as in the model of Park &z Lee! (1997, see also Yi et al 



19991 ) . Our composite models with helium-enhanced subpopulations can reproduce the UV- upturn 
at the mean age of 11 Gyrs, in the age scale where the inner halo GCs of the Milky Way i s 12 Gyrs 
old. If the age of the oldest GCs in the Milky Way is more like 13 Gyrs ([Potter et al.ll2010i ). our age 
for the giant elliptical galaxies would also increase by 1 Gyr. These ages for giant elliptical galaxies 



1 For the dwarf elliptical galaxy like M32, where the observed UV upturn is weak (FUV — V = 7.3 in AB 
magnitude system), all of the HB population is predicted to be cooler than 20,000 K. Brown et al. 1 200oh found from 
their HST/STIS photometry of M32 that stars passing through the HB at T c g > 8,500 K comprise only a small 
fraction (about 7%) of the total HB population. This fraction is roughly consistent with that (~ 5.5%) predicted 
from our models for M32. 

2 Most of the helium-rich HB stars in our models are too hot to have any significant effect on metal lines, such as 
Ca II H + He/Ca II K index. The Ca II H + He/Ca II K index of our model f or giant elliptical galaxy (presented in 



Figure 3) is 1.15, which is consistent with the observation (~ 1.16; Rose 1985) to within the error (± 0.03) 
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(11 - 12 Gyrs) are consistent with the value e stimated from the Balmer line indices to within the 
error (jlYager et al.1 hood : iThomas et al.ll2005l : llYaeer et alibood : braves et al.lbood 'l. 



In Figure 4, we have plotted the Mg b line strength against FUV — V color for the sample of 



quiescent ETGs from lBureau et al.1 (120111 ). Superposed grids are our composite models with (blue) 

The observed ETGs show clear correlation 



Burstein et al.l ([1988) relation. Our composite 



and without (red) helium-enhanced subpopulations 
between FUV — V and Mg b that is analogous to the 
models predict that the strength of UV-upturn is controlled by the fraction of helium-enhanced 
subpopulation, the mean age, and the mean metallicity of underlying stellar population. Among 
these three factors, the fraction of helium-enhanced subpopulation appears to be the most effective 
factor determining the strength of UV-upturn. For example, our 11 Gyr models with and without 
helium-enhanced subpopulations would explain most of the color spread in FUV — V. If all ETGs 
in this sample, however, contain similar fraction of helium-enhanced subpopulations, an age spread 
spanning ~ 6 Gyrs would be required to reproduce the observed spread in FUV — V color, in the 
sense that bluer galaxies are older. This is because the mean temperature of HB stars decreases 



as the age of stellar population gets younger (see iLee et al 



1994 



spread of this magn i tude is also indicat e d from the Balmer line dating of ETGs (jTrager et al. 



Thomas et al.ll2005l ; iTrager et al. 



2008 



Park & Led 


1991). An age 


of ETGs ( 


Tracer et al. 


2000; 



Graves et al.ll2009l ). which suggests that some age spread 



is indeed responsible for the spread in FUV — V color. The obs erved fading of FUV flux with 
look-back time for the bright cluster elliptical galaxies in z < 0.3 (|Ree et al.l 120071 ) also supports 
this possibility. Note from Figure 4 that our models both with and without helium-enhanced 
subpopulations predict fading of FUV flux with increasing look-back time. 



4. DISCUSSION 

We have demonstrated that the presence of relatively metal-poor and helium-enhanced subpop- 
ulations in ETGs can naturally reproduce the observed UV-upturn phenomenon, without invoking 
unrealistically old ages. Although the presence of helium-enriched stars appears to solve many of 
the problems associated with the UV-upturn phenomenon, it is important to note that the origin 
of this helium enhancement is not fully understood yet. Consequently, it is not clear yet whether 
the proposed mechanisms to produce helium-enhanced stars in GCs will also be important in the 
galactic scales. Nevertheless, several lines of evidence do suggest that a minority population of 
helium-enhanced stars would be present also in ETGs. First, the two well defined remaining local 
building blocks in the Milky Way, the u) Cen and the Sagittarius dwarf galaxy (includi ng M54), 



are all characterized by very extended HB with helium-enhanced extreme blue HB stars (jLee et al 



3 Our models were constructed with [a/Fe] = 0.3, while some galaxies in Figure 4 have lower values of [a/Fe]. 
If [a/Fe] = 0.1, the FUV — V color gets redder by small amount (~ 0.1), while the Mg b index is more affected 
(decreases by ~ 0.4 A). Therefore, the model grids in Figure 4 shift mostly downwards as [a/Fe] decreases, and this 
should have only little effect on the predicted age spread among sample galaxies. 
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2005; 


Piotto et al. 


2005; 


Bellini et al. 


2010; 


Sieeel et al. 


2007: 



2011, in prep.). Second, the orbital kinematics of Milky Way GCs with extended HB, which is 
indicating the presence of helium-enhanced subpopulation, are distinct from normal GCs, and are 
fully consistent with the hypothesis that they are the remaining relics of the early bu i lding blocks 



predi cted in the hierarchical merging paradigm of galaxy formation (|Lee et al.l 120071 ; iBekki et al 
20071 ). Finally, as discussed above, when the observed trend of FUV — V color with metallicity 



for the UV-strong GCs (presumably with the extended HB) in the giant elliptical galaxy M87 is 
combined with the expected MDF of a model giant elliptical galaxy, the FUV — V color similar to 
that observed in giant elliptical galaxies is reproduced. This is again consistent with the building 
block origin for these GCs, according to which the helium-enhanced subpopulations in ETGs would 
have been supplied by these early galaxy building blocks. 

If normal ETGs are prevailed by helium-enhanced populations, as suggested in this Letter, they 
would have impacts not only on UV, but also o n other waveleng th regime affected by relatively hot 
blue HB stars, such as B a lmer absorption lines (|Lee et al.ll2000l ). In this respect, it is important to 
note that lSchiavon et al.1 (|2006l ) found that their models (without helium-enhanced subpopulation) 
of passive evolution can not reproduce the unexpectedly strong Balmer line indices of ETGs at z ~ 
0.1 and z ~ 0.9. It would be interesting to see how our models with helium-enhanced subpopulations 
are matched with these observatio ns. The well kno wn discrepancy between the H/3 and higher order 
Balmer lines ("Balmer mismatch" ; ISchiavonll2007l ) is also an important issue, for which the helium- 
enhanced populations could play a role. Our forthcoming paper (Chung et al. 2011, in prep.) will 
discuss these issues in more detail, together with their connection with the UV-upturn phenomenon. 
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Fig. 1. — The synthetic HR diagrams and corresponding spectral energy distributions for simple 
stellar populations. Upper panels are models for normal helium abundance (Y = 0.23), while 
bottom panels are for helium-enhanced (Y = 0.33) population. The abbreviations MS, RGB, HB, 
and PAGB refer the flux contributions from main-sequence, red giant branch, horizontal branch 
(including post-HB progeny), and post asymptotic giant branch stars, respectively. 
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Fig. 2. — (a) Comparison o f our model s (soli d lines) with GCs in M87 (circles) and t he M ilky Way 
(triangles). The data are from ISohn et al.1 (|2006) for GCs in M87 and from iDorman et al.1 (1995) for GCs 
in the Milky Way, and the FUV — V color is in the AB magnitude system. Bared arrows are upper limits 
for the metal-rich GCs in M87. The red line is a model with only normal helium population (Y = 0.23; 
t = 12 Gyr), while the blue line is a model when the helium-enhanced subpopulations (Y = 0.33 + 0.38; 
t = 11 Gyr) are added to the normal helium population (Y = 0.23; £=11 Gyr) using the fraction in 
Figure 2b. Some well-known Milky Way GCs with extended HB (w Cen, M54, NGC 6388, and NGC 6441), 
together with 47 Tuc, are labeled. [Fe/H] values for ui Cen and M54 are weighted mean values of multiple 
populations in these GCs. (b) The fraction of helium-enhanced population calibrated to reproduce FUV — V 
colors of UV-strong GCs in M87 at given metallicity. In our models (the blue line in Figure 2a), helium- 
enhanced subpopulati ons are added to the norm al helium populations with this relation. ( c) The metallicity 
distribution function (jKodama &: Arimotdll997l ) adopted in our composite model for giant elliptical galaxy 
and the contribution from helium-enhanced subpopulations. Considering the metallicity spread of composite 
model peaked at [Fe/H] = 0.3, the actual contribution from helium-enhanced subpopulation is only ~ 11%. 
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Fig. 3. — The comparis on of observed SEDs of NGC 4552 and NGC 4649 (data from lBurstein et al 
19881 and lYi et al.lll998l ) with our composite models. The composite models presented here are based 
on the series of SSP models such as those displayed in Figure 1, with the metallicity distribution 
function in Figure 2c. Our models with helium-enhanced subpopulation can reproduce the observed 
far-UV upturns at the age of 11 Gyrs. 
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Fig. 4. — The (FUV — V)q colo r vs. Mg b correla tion for the sample of quiescent early-type 
galaxies having H/3 < 1.9 A from bureau et al.1 ( 201ll ). The blue and red lines are our composite 
models with and without helium-enhanced subpopulation, respectively, where the iso-metal lines 
are connected by dotted lines. The values of (Z) and age (Gyr) are indicated. 
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Table 1. INPUT PARAMETERS ADOPTED IN OUR COMPOSITE STELLAR 
POPULATION MODEL FOR NGC 4552 



Parameters Helium-enhanced population Normal Helium population 

Initial mass function Salpeter (x = 1.35) Salpeter (x = 1.35) 

a-elements enhancement, [a/Fe] 0.3 0.3 

HB mass dispersion, cm 0.015 0.015 

Reimers' mass-loss parameter, r) 0.63 0.63 

Absolute age, t (Gyr) 11.0 11.0 

Mean metallicity, {[Fe/H]) -0.76 -0.02 

Helium abundance, Y 0.33 + 0.38 0.23 

Population ratio (%) 11 89 



